ZEROES OF THE SPECTRAL DENSITY OF THE 
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Abstract. We consider the Schrodinger operator L a on the half-line with a periodic 
background potential and the Wigner-von Neumann potential of Coulomb type: 
£SS^g±i) _ \i j s known that the continuous spectrum of the operator C a has the same 
band-gap structure as the free periodic operator, whereas in each band of the absolutely 
continuous spectrum there exist two points (so-called critical or resonance) where the op- 
erator C a has a subordinate solution, which can be cither an eigenvalue or a "half-bound" 
state. The phenomenon of an embedded eigenvalue is unstable under the change of the 
' boundary condition as well as under the local change of the potential, in other words, 

it is not generic. We prove that in the general case the spectral density of the operator 
C a has power- like zeroes at critical points (i.e., the absolutely continuous spectrum has 
pseudogaps). This phenomenon is stable in the above-mentioned sense. 
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1. Introduction 



It is well-known that the Schrodinger operator on the half-line with a summable potential 
has purely absolutely continuous spectrum on R + [28]. It is also known [12J that if the 
potential is square summable, then the absolutely continuous spectrum still covers the 
positive half-line, see also the work [T7]. In this case however there is no guarantee that 
the positive singular spectrum is empty. Moreover, there are explicit examples of potentials 
f25[ I2T], which produce dense point spectrum on R + . In a sense, the simplest example of 
■ the potential producing positive eigenvalues is the so-called Wigner-von Neumann potential 
csm(2uix) ^ se j£_ acij oi nt operator £ ,« given by the differential expression / : = — ^ + 
esin(2oja:) &n( j boundary condition ^(0) cos a — ip'(0) sin a = 0, acting in L 2 (IR + ) on the 
domain 



Dom£ , a = jv> e L 2 (R + ) n Hf oc (R + ) : l ip e L 2 (M+),^(0) cos a - ^'(0) sin a = oj 

may have an eigenvalue at the point u 2 [30J, while the rest of M + is covered by the purely 
absolutely continuous spectrum [2]. Due to the subordinacy theory [13], this can be seen 
from the behavior of generalized eigenvectors (i.e., the solutions of the spectral equation 
l ip = Xijj): for A G M+\{w 2 } there is a base of them with asymptotics e ±l ^ x + o(l) as 
x — > +oo. For A = bj 2 the base is different: '(sin(a;a;) + o(l)) and x~^ (cos(ux) + o(l)). 
This point is called the resonance (or critical) point due to the change of the type of 
asymptotics of generalized eigenvectors. One of the solutions (up to a constant) of the 
equation l ip = u 2 ip is subordinate. If on top of that it belongs to L 2 (M + ) and satisfies 
the boundary condition at the origin, then it represents an eigenfunction of £o,a- It is 
clear that the effect of appearance of positive eigenvalue is highly unstable. It happens 
with probability one in a suitable sense: the eigenvalue disappears if one changes slightly 
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the boundary condition or adds a summable (or even compactly supported) perturbation 
to the potential. If the notion of resonances makes sense for the operator (i.e., if the 
kernel of the resolvent admits the analytic continuation through the continuous spectrum), 
this eigenvalue may become a resonance under such perturbation. The set of summable 
functions which one could utilize as perturbations of the potential preserving the eigenvalue 
at the point u 2 was studied in [TO]. On the other hand, the phenomenon of the change of 
the type of asymptotics of generalized eigenvectors at a point of the absolutely continuous 
spectrum is stable under summable perturbations of the potential and does not depend on 
the boundary condition. It is therefore meaningful to consider objects which are a more 
or less stable in this sense: one can study the Weyl-Titchmarsh function m or the spectral 
density p' (the derivative of the spectral function of the operator), which are related by 
the equality 

p'(A) = - Imm(A + iO) for a.a. Ael. 

7T 

The behavior of the Weyl function near resonance points has been studied by Hinton-Klaus- 
Shaw [151 118] and Behncke [21 El H] for Schrodinger and Dirac operators with potentials of, 
in particular, Wigner-von Neumann type. 

In the present paper we consider the differential Schrodinger operator C a with the po- 
tential that is the sum of the following three parts: a periodic background q, a potential 
of Wigner-von Neumann type csm ( 2 ^ E+<5 ) and a summable part q\. The operator is defined 
by the differential expression 

. d 2 . . csm(2ux + 5) . . 

(1) t;= -_ + gW+ ^ + i > +qi[x) 

and the boundary condition ip(0) cos a — ip'(0) sin a = 0, a G [0; it). It acts in L 2 (M + ) on 
a suitable domain. We assume that the function q has the period a and is summable over 
this period: q 6 -£q(0; a). The parameters c and 5 are real constants. The operator C a is 
then self-adjoint in L 2 (M + ). 

The geometry of the spectrum in the case of periodic background differs from the case 
q(x) = 0. It is well-known (see, e.g., [22]) that the absolutely continuous spectrum of the 
operator given by the expression I on the whole real line coincides with the spectrum of 
the corresponding periodic operator on the whole real line, 

d 2 

(2) C per = -— +q{x), 

i.e., has a band-gap structure: 

oo 

v(£per) ='■ [_J([^2j; A*2j] U \l 1 2j+l] ^2j+l] ) ) 

where A < Po < p\ < Xi < A 2 < pi < pz < A 3 < A 4 < ... It is clear that the absolutely 
continuous spectrum of C a coincides set- wise with a(C per ), although a(C per ) has multi- 
plicity two, whereas a(C a ) is simple; note that C a is considered on the half-line. In every 
band there exist two resonance points z^- i+ and The type of asymptotics of gener- 

alized eigenvectors at these points is different from that in other points of the absolutely 
continuous spectrum. Subordinate solution appears and thus each of the resonance points 
can be an eigenvalue of the operator C a as long as this solution satisfies the boundary 
condition and belongs to L 2 (IR+). The exact locations of the points Vj t ± are determined by 
the "quantization conditions" [22J 

k(u ji+ ) = TT (j + 1 - {^}) , k(u jt .) = TT (j + |^}) , j > 0, 
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where k(X) is the quasi-momentum of the periodic operator C per and {•} is the standard 
fractional part function. The above-mentioned condition w ^ ff guarantees that resonance 
points firstly do not coincide with the band boundaries and secondly that they do not glue 
up together. 

The main goal of the present paper is to analyze the interplay between the periodic 
structure and the "singular" Wigner-von Neumann perturbation. The main result is The- 
orem [2] of Section [51 which essentially tells us that the spectral density of the operator C a 
has power-like zeroes at each of the resonance points. Let ip + (x,\) and ip-(x, A) be the 
Bloch solutions of the periodic equation —i])"(x) + q{x)^(x) = \ip(x). Let (f a (x, A) be the 
solution of the Cauchy problem 



(3) l<p a = \ip a , <f a (0) = sin a, <p' a (0) 



cos a. 



Denote by W{i/) + (\),ip-(\)} = W{ip + (-,\),ip-(-,\)} the Wronskian of two Bloch solu- 
tions. 

Theorem. Let qi £ Li(M + ) ; u ^ and p' a (X) be the spectral density of the operator C a . 
Let the index j > 0. If a is such that the solution (p a (x, Vj,+) of ([3]) is not a subordinate 
one, then there exist two non-zero limits 

Km ^ 



\-^v, +±o M 



Analogously, if a is such that tp a {x,v^J) is not subordinate, then there exist two non-zero 
limits 



lim 



gjcj 



\X-v jr 



a\W{-4, + {u^ + )rt_{»j,+)}\ 



f tpl(t,v :j -)e2<«> t dt 



(i 



/ i) 2 ± (t,v cr )e 2iult dt 
o 

ows from |22| anc 



Note that the condition of non-subordinacy of (p a is indeed satisfied in the generic case. 
Consider one of the critical points v cr . There is a unique a £ [0; 7r) for which ip a is 

subordinate, denote it by a cr . Denote also B := ,„, f , , H , , — viT 

where the choice of sign coincides with the one in v CT = Vj t ±. It fol 
the analysis below that ip acr (x, v CT ) = x~^{c + ip + {x, z/ cr ) + c-ip-(x, v cr ) + o(l)) as x — > +oo 
where ip+ii^- are the Bloch solutions of the periodic equation —ip"(x) + q(x)ip(x) = \ip(x) 
and c + ,c_ are some constants. Therefore v cr is an eigenvalue of C acr iff 8 > \. At the 
same time, for 8 > \ and a ^ a cr the order of the zero of p' a (X) at the point v CT is greater 
than one. This fact is in exact correspondence with the result of Aronszajn-Donoghue PQ, 
by which if the point v is an eigenvalue of the operator C acr , then the integral J R r£z~^ is 
finite for every a ^ a cr . 

It should be mentioned, that this result in the particular case of zero background pe- 
riodic potential q follows from the work of Hinton-Klaus-Shaw |15j . They considered the 
differential Schrodinger operator on the half-line with the potential, which is an infinite 
sum of terms of Wigner-von Neumann type plus a rapidly decreasing term analogous to our 
qi but decaying faster. They studied the behavior of the Weyl-Titchmarsh function near 
the critical points in the case a = a cr and 8 < \ , i.e., when the solution (p a is subordinate 
but is not an eigenfunction (the so-called half-bound state). The corresponding result for 
the case a ^ a cr and any 8 > was a by-product of their analysis. Nevertheless, the most 
essential part of the problem in their work is similar to the one in our case. However, our 
approach is different from that of papers [151 dH] , see the detailed discussion in Section El 
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Our goal was to elaborate a more general approach (reduction to a model problem), which 
would enable us to apply it without any significant changes to different operators. 

Zeroes of the spectral density of the Schrodinger operator divide the absolutely continu- 
ous spectrum of the operator into independent parts. This phenomenon is called a pseudo- 
gap and has a clear physical meaning. Eigenvalues embedded into the continuous spectrum 
have been observed in experiment [H] . Operators with Wigner-von Neumann potentials at- 
tracted attention of many other authors, e.g., El HS1 HSl SI El El EH 1201 1221 - 

The paper is organized as follows. In Section [21 we state the Weyl-Titchmarsh type 
formula for the operator £ a which was proved in [23J and serves as our main tool in the 
work with the spectral density. In Section [31 we write the spectral equation for the operator 
C a in an equivalent form of a discrete linear system and restate the Weyl-Titchmarsh type 
formula in terms of the asymptotic behavior of solutions of that system. In Section HI 
we study this discrete system (which can be regarded as a model problem). This can be 
viewed as an independent task. In Section \5\ we use the results of the preceding analysis 
to find the asymptotics of the spectral density of the operator C a near the critical points. 
In Section El we make final comments characterizing our method in comparison to the one 
utilized in the work [15] and give a discussion of a few examples of its applicability. 



2. Preliminaries 

Spectral properties of second-order differential Schrodinger operators are related to the 
asymptotic behavior of their generalized eigenvectors |13j . One of illustrations of this 
relation is the Weyl-Titchmarsh type formula for the spectral density in terms of asymptotic 
coefficients of the solution to the spectral equation that satisfies the boundary condition. 
The spectral equation lip = Xip can be considered as a perturbation of the periodic equation 
—ip"(x) + q(x)ip(x) = X%p{x). 

Let us denote by d := {Xj,fij, j > 0} the "generalized boundary" of the spectrum of 
the periodic operator C per (it can be not exactly the boundary of the set a(C per ), because 
the endpoints of the neighboring bands can coincide). We will use the following result 
(Weyl-Titchmarsh type formula). 

Proposition 1 ( [23j ) . Let ^ ^ ff o,nd q\ G Li(R+). Then for every fixed 
X G o~(£ per )\(d U {vj i+ ,i/j^,j > 0}) there exists a non-zero constant A a (X) depending 
on X such that 



, . (p a {x, X) = A a (X)ip-(x, X) + A a (X)ip + (x, X) + o(l) as x -)■ +oo, 

1 ' ip'Jjc, X) = A a (X)ip'_(x, X) + A a (X)ip' + (x, X) + o(l) as x -> +oo, 

and the following equality holds: 
(5) P'M 



2tt\W{iP + (X),iP_(X)}\ K(A)P' 

where p' a {X) is the spectral density of the operator C a . 

This result is a generalization of the classical Weyl-Titchmarsh (or Kodaira) formula 
[29l [19] . A variant of this formula for the Schrodinger operator with Wigner-von Neumann 
potential without the periodic background (q(x) = 0) follows from the results of [2^1 16]. In 
the case of discrete Schrodinger operator with Wigner-von Neumann potential an analogous 
formula is also known, see [HI [To] . 
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3. Discretization 

In this section we pass over from the spectral equation lip = Xip to a specially constructed 
discrete linear system essentially equivalent to lip = \ip and rewrite the Weyl-Titchmarsh 
type formula expressing the spectral density in terms of solutions of this system. 

Step 1. Variation of parameters and discretization: perturbation of the mon- 
odromy matrix. First we perform a transformation of the spectral equation to the differ- 
ential system in order to "get rid" of the periodic background potential (in fact, this simply 
is a variation of parameters). Consider A G a(C per )\d and equation lip = Xip. Define new 
vector- valued function r\ [x) by the following equality: 

*f>( x ) \ _ ( ^-(z,A) Mx,X) 
iP'{x) ) \iP'_{x,\) iP' + (x,\) 

By a straightforward substitution one has: 

(6) r]'(x) = L(x,X)r](x), 

where 

L{Xt A) ._ csH Ti +8) + ( -m%, \)M*> ^ -r + (x, x) 



rj(x) 



W{ip+(X),ip4X)} V CO,A) ip+{x, X)ip-{x, X) 

Denote by $(zq, x, X) the fundamental matrix for the system (i.e., the matrix solution of 
(j§D satisfying ^(x , x , X) = I). Denote also by M n (X) the monodromy matrix correspond- 
ing to the shift by the period a of the background potential, M n (X) := $(a(n — 1), an, A). 
Instead of the solution t](x) of the differential system ([6]) consider the following sequence 
{t« n }™ j of vectors from C 2 : 

w n := i](a(n - 1)), 

which obviously solves the discrete linear system w n+ \ = M n (X)w n . In what follows, we 
deal with this system transforming it to a simpler form. Using the standard perturbation 
methods (see for example [22]) we can determine the matrix M n (X) up to a summable (in 
n) term. 

We have passed from the continuous variable x to the discrete variable n. Later, in 
Section H] we will return to continuous variables considering a new one, y. There is no 
strict necessity in such a trick. One could think that the discretization determined by the 
period a of the background potential q helps to get rid of this potential. However, this is 
not completely true: this role is taken by the above variation of parameters transformation. 
One could do simpler: the coefficient matrix of the differential system (JS]) is the sum of 
two terms. The first term has a factor ^j-j- multiplied by an infinite sum of the exponential 
terms coming from the Fourier decomposition of the periodic parts of Bloch solutions. 
Each of these terms can "resonate" (become constant in x) for certain values of A. This 
happens exactly at the resonance points Vj,±- The second term is a summable matrix- valued 
function. Non-resonating exponential terms can be eliminated using the uniform Harris- 
Lutz transformation in a fashion of [23J. This approach would lead to a simple differential 
model system equivalent to the equation lip = Xip. However, we use the approach of the 
discretization, which brings some extra technical difficulties and complicates the notation. 
The reason is that we obtain as a result a discrete model system, which can serve wider 
needs. In particular, it is possible to consider the discrete Schrodinger operator with 
the discrete Wigner-von Neumann potential which has two critical points on the interval 
[—2; 2] (covered by the absolutely coninuous spectrum of this operator) and to prove that 
its spectral density has zeroes of the power type at the critical points. See Discussion for 
the exact formulation of this result. We plan to prove it in a forthcoming paper. 
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Let us introduce the following notation for uniformly summable sequences. Let Rn(X) be 
2x2 matrices for n G N and A G S (where S is an arbitrary set). We write {R n {X)}^ =1 G 
^(S), if there exists a sequence of positive numbers {r n }^ =1 G I 1 such that for every A G S 
and n G N, ||i? n (A)|| < r n . We start with a simple asymptotic formula for the monodromy 
matrix: 

Lemma 1. For every A G a(C per )\d, 
M m - — r" csin(2a;t + £) 



-V + (t,A)V-(t,A) -^(t,A) 
^l(t,A) V+(*,A)V-(t,A) 



where {R^\X)}^ =1 G Z x (i^) /or every compact set K C a(C per )\d. 

Note that it follows that for every n G N the matrix i?^(A) depends continuously on 
A G a(C per )\d. 

Proof. We give only a sketch of the proof, because it is rather standard. Fundamental 
matrix for ([6]) satisfies the integral equation 

(7) $(x ,x, A) = 1+ / L(t,X)$(x ,t,X)dt. 



$Oo,x,A) = 1+ / L(t,A)$(x ,t, 



To study M n (A) put x = a(n— 1) and consider Volterra integral operator V n (A) in Banach 
space C([0; a], M 2,2 ) (of 2x2 matrix functions, with any matrix norm) defined by the rule 

V n (A) : u(x) t-> / L(a(n -1) +t,X)u(t)dt. 
Jo 

Using Neumann series and direct estimate of the norm we can show the existence of (J — 
V„) _1 and the following estimate for its norm: 



an 



(8) ||(/ -y^A))- 1 !! <exp / \\L(t,X)\\dt 

\Ja(n-l) j 

Returning to (J7J) we can write the matrix equality 

$(a(n - 1), a(n - 1) + x, A) = ((/ - V n (A))- 1 /)(x) = 

= / + (V„(A)/)(x) + {Vl{X){I - V n {X))- x I){x) 

and so putting x = a we have: 

pan 

M n (A) = / + / L(t, X)dt + (V 2 n (X)(I - V„(A))~ 1 /)(a) . 

Ja(n-l) * v ' 



i(n-l) 

= :i?l 2) (A) 



The remainder i?„ (■) is continuous in a{C per )\d for every n and satisfies the uniform 
estimate 

rV\X) = 0\(J^ \\L{t,X)\\dtj\ as b ^ oo. 

Fix the compact set K C o~(£ per )\d. One can always choose Bloch solutions so that their 
Wronskian is analytic functions and does not vanish on o~(£ per ). From the properties of 
Bloch solutions it follows that there exists C\(K) such that for every A G K and x, 

1 ( -if)+{x, A)V>-(x, A) - ipl(x, A) 

1U{^+(A),V-(A)} V ^+(x,X)i;.(x,X) 



< Cl (K). 
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A straightforward estimate yields: 

r \\L{t,\)\\dt< Cl (K)(\c\\n( an+l )+r \ gi (t)\dt), 
Ja(n-i) V \an + l-aj J a(n ^ 1} J 

and thus {i?i 2) (A)}~ =1 G l\K). Now consider the expression for §™(n-i) L(t, X)dt. Fix the 
parameter t in the denominator, putting it equal to an (this gives precisely the expression 
for the integral in the assertion of the lemma). The difference is: 
(9) 

csm(2cot + 5) ( 1 1 \( -<0 + (t,A)V-(t,A) -ipl(t,\) \ df 



which admits an estimate by ^(^ONfo+i) ^he total error term i?^ is the sum of Rn^ and 

J n(an+l—a) 

the expression and therefore it possesses the required properties. This completes the 
proof. □ 

Step 2. Fourier decomposition of Bloch solutions. Rewrite the expression from the 
second term of the formula for M n (A) given by Lemma [TJ 

If csin(2u;t + 5) ( -ip + (t, \)i^-(t, A) A) \ , 

1 J «Vi) W{iMA),lMA)} V ^M) ^+(*,A)V»-(t,A) j 

using Fourier decompositions for periodic parts of Bloch solutions. Let {&j l "(A)}~_ DO , 
{&r(^)}£-oo anc ^ {^(^)}£-oo be Fourier coefficients defined by identities 

+oo +oo 

(11) 4 + (x,\)1,-(x,\) = £ W)*"* 1 *, V£(*,A) = £ 6f (A)^^^. 

These Fourier coefficients and their derivatives with respect to A are locally uniformly 
bounded in a(C per )\d and obey the locally uniform (on the same set) estimate 0{l/l 2 ) as 
/ — > oo. This fact is rather standard, see, e.g., [23] for the details. Since ip + and ■?/>_ are 
complex conjugate on a(C per )\d, one has: 



6+ (A) = 6I,(A), 6j(A) = 6_ { (A) for every I £ Z and A G a(£ per )\d. 

Substituting identities ffTTj) into the expression ffTOj) and changing the order of summation 
and integration (which is possible due to the properties of Fourier coefficients mentioned 
above), one obtains the following result. Expression ( 1TU1) equals to 

#(A) M ad) (A) 



^ ad) (A) -#(A) 



, and the monodromy matrix can be written in the form 



1 / A d '(\) ftf*(A) 



(12) M„(A) = / + ± ^ +fl<f(A). 



77 



/^ ad) (A) -M d) (A) 



Here we have introduced the following notations: 



a (\\ ._ ce i(S-au) „ , . sin(aa;) 

POl^J •— 2W{V+(A),V-(A)} ^ Ul \ A ) irl+aui ' 

(13) /= + ~^ 

o /\\ ce i(6-(k(X)±au)) , + , . , S m(k(\)±au) 

l=—oo 

and 

._ a i \ 



( 14 ) M 



/3r(A):=/3 (A)e 2 ^-/3o(A)e' 
(3 { n d \\) := / 3 + (A) e 2j ( A; ( A ) +a ^ n + f3_(\) e MkW-au)n_ 
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Step 3. Elimination of non-resonant terms by Harris-Lutz method. For every 
fixed A G o-(Cper)\d the entries of the matrix I ^zJ ~ ^ ^ n ,1 ^ \ are linear eonibma.- 



M ad) (A) -/?r(A) 

tions of four exponential terms: e 2iawn , e " 2iaum , e 2i{k{ ^ +a ^ n and e 2i(k(A)-ow)n The first two 
of them do not depend on A and since oj ^ |^ they do oscillate. The third and fourth terms 
do depend on A. They oscillate, if k(X) +aw ^ 7rZ and k(X) - aw ^ 7rZ, respectively, other- 
wise they are constant (we call resonance this matching of the quasi-momentum k(X) and 
the frequency u>). Oscillating exponential terms can be dropped with the help of Harris- 
Lutz transformation and do not affect the asymptotical type of solutions of the system 
w n+ i = M n (X)w n . Therefore (see the rigorous proof of this fact below) if A is not one of 
the resonance points Vj,±, which are defined by conditions 

fcfe) := - (j + 1 " { f }) , H"j,-) := -r (j + {^}) ,j>0, 

then every solution of the system w n+ \ = M n (X)w n has a limit as n — > oo. Note that due 
to the condition to for every j > resonance points z/j >+ 7^ Z/, • _ do not coincide with 

the endpoints of the j-th spectral band. 

Pick a spectral band and one of the two critical points in this band, e.g., for some 
index j. To simplify the notation, we write u cr : = z/ Ji+ . Consider an open neighborhood 
U cr of the point v cr such that it contains neither the critical point Vj_ nor the endpoints 
of the band fij and Xj. In what follows we assume that v cr and U cr are fixed and we drop 
the indices j, + in most cases. To remove the oscillating terms from the system we need 
the following elementary lemma belonging to the domain of the mathematical folklore. 

Lemma 2. For every real £ ^ 2nZ and n G N, 



E 



m 



< 



n sin 



Proof. 



„irrt£ 



m 



E 



m 



m + 1 



71 



^Ef 1 - 

rn — n ^ 



m + 1 



12 

71 71 



□ 



Since |e^ — 1| =2 |sin ||, this argument completes the proof. 

By a Harris-Lutz transformation uniform in U cr it is possible to eliminate non-res- 
onating exponential terms and to stabilize coefficients at resonating terms (i.e., make them 
independent of A). Doing this we can provide for the uniform summability of the remainder. 
We formulate this argument in the following lemma. 

Lemma 3. There exists a sequence of matrices {Qn\X)}^ =1 defined in U cr such that 
Qn\X) = O (-) as n — >■ 00 uniformly in U cr and 



(15) exp(-Q« 1 (A))M„(A) exp(QW(A)) 



n 





-2i(k(\)+au>)n 



a l \ 2i(k{X)+au))n \ 



w^/i some {Rn\X)}™ =l G / 1 (f/ cr ) snc/i £/ia£ .Ri 3 ^(-) continuous in U cr for every n. 
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Proof. We follow the scheme of [5] and need to take care of the uniformity only. The 
explicit formula for Qn ^ is 

oo _ 

(16) Q« (A) :=-£)- 



- A d) (A) 



Lemma [2] yields immediately the estimates 



E 



0m (X) 



m 



< 



2|/3o(A)| 



n sin aw 



and 



E 



A° J) (A)-^ + (^.)e K ( i < A >+""'> 



< |/3 + (A)-/3 + (z/ cr )| | 



l/MA)| 



n| sin(/c(A) + au;)| n| sin(fc(A) — aw) 



Functions (3q and /3_ are continuous and the denominators sin aw and sin(fc(A) — au) are 
separated from zero in U cr . On the other hand, the denominator sin(fc(A) + au) has the 
only zero at the point v cr , which is simple and thus is compensated by the zero of the 
numerator. Therefore the estimate Qn\X) = O (-) holds and is uniform in U cr . Using the 
obvious property 

& d \\) P { n ad \X) - /3 + (z/ cr )e 2l « A )+^ 

& ad \\) - f3 + (u cr )e- 2i ^ +a ^ n 



QSi(A)-gi 1) (A) 
one obtains: 



1 

n 



M n (A) 



1 



(3 + (u cr )e 2i ^ +aw ^ n 

+ Qi 1 !i(A)-Qi 1) (A) + 4 1) (A). 

Multiplying by exp(Qn\\)) from the right; by exp(— Qn+i(X)) form the left, expanding 
exponents and absorbing the terms of the order 1/n 2 in the remainder, we have: 



exp(-Qi 1 ) 1 (A))M„(A)exp(g^(A) 



>(!)/ 



n V (3 + {u cr )e- 2i ^ +a ^ 



=«1 3) (A) 



where the estimate O (-^) is uniform in U cr . It is clear that R n 3 \X) is continuous in U cr 
for every n. This completes the proof. □ 



Step 4. Reduction to the model problem. After the Harris-Lutz transformation 

{Wn}n=l ^ {Wn\n=l with 

w n := exp(— Qn \X))w n the system w n +\ = M n (X)w n is reduced to the system 

/3 + (u cr )e 2i ^ +a ^ n 




n V, (3 + (is cr )e- 2i(k ^ +aw ^ n 
At the critical point it takes the following form: 



+ i2?>(A) 



w„. 



Wn+l 



1 







/3+(^cr) 

n V /3+(z/ cr ) 
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gfarg P+(v cr ) 
g-f arg£ + (i/ cr ) 



P+M 

(3 + {u cr ) 



e f arg/3+(i/ cr ) ^gf arg/3+(z/ cr ) 
g-f arg/3 + (V cr ) _^g- 5 arg /3 + (^ cr ) 



1 
-1 



Consider the following sequence {v n }^ =l instead of {-uv}^: 



gf arg/3+(i/ cr ) ^ e |arg/3 + O cr ) 
g-f arg/3+(z/ cr ) _jg-f arg/3+(i/ cr ) 



Due to ffl"5l) the system to n +i = M n (X)w n is equivalent to: 
(17) 



■v, 



n+l 



T , P_ ( cos(2(k(X) + au)n) sin(2(/c(A) + aw)n) \ d(4)/\\ 
+ n I sin(2(fc(A) + acu)n) - cos(2(fc(A) + aw)n) ' " r " ( ' 



where the sequence {k£\\)}™ =1 G l l (U cr ), the functions -Rn 4-I (-) are continuous in U cr and 

k'l 



>(4), 



2a|iy{^ + (z/ cr ),V_(i/ cr )}| 



Replace the parameter A on the set U cr by the new small parameter 
(18) e := 2(A;(A) - k{v cr )) = 2(k(X) + aw) - 2tt (j + 1 + 



aw 



where |_'J is the standard floor function. The set of values taken by e is U := {2{k{\) — 
fc(z/ cr )),A G U cr }. By the property of the quasi-momentum (that k(X) G [7rj;7r(j + 1)] in 
j-th spectral band) and the condition to ^ which guarantees that the critical point is 

in the interior of the spectral band, we have U C (— 2tt; 2tx). Denote R n (e) = B£ (A) for A 
corresponding to e according to (fi~8l . System (TTTj) then reads: 



(19) 



1 + 



/3 



n 



cos (en j 
sin (en) 



sm{en) 
- cos(en) 



+ Rn(e) 



The aim of this section is to rewrite the Weyl-Titchmarsh type formula in terms of the 
solutions of the system (TT9l . Proposition [1] deals with the solution ip a of the spectral 
equation for the operator C a . Combining all the transformations described above and 
denoting the result by S : ip(x) h-> {vn}^^, one ends up with the following model image 
of the initial solution ip a : 



(20) v a , n (e) = v a , n (2{k{\) - k(vcr))) 



g§ arg/3 + (i/ cr ) arg/3 + (u cr ) 
e ~5 arg/3 + (i/ cr ) _ ? 'g-| arg£ + (v cr ) 
-1 



ip_(a(n — 1), A) V ; +( a ( ri ~ I); A) 



xexp( (A)) l ^( Q ( n _ 1 ) )A ) ^( a ( n _i),A) 
This is obviously a solution of the system ( TT9]) and it is continuous in {7 for every n. 



<p a (a(n - 1), A) 
^ a (a(n-l),A) 



Lemma 4. For ever?/ s G t/\{0} i/iere exists a limit lim u a n (e) 7^ 0, which is continuous 
in e G C/\{0} as a vector-valued function of e. The spectral density of C a equals 

2 



(21) /4(A) 



tt|W{V + (A)^-(A)}| lim u Qin (2(fc(A)-A:(i/ cr ))) 



a.e. in C/ r , 



ZEROES OF THE SPECTRAL DENSITY 11 

Proof. A straightforward substitution of the transformation H yields this result. Indeed, 
asymptotics of the solution ip and of its derivative @ together with the boundedness of 

i[>-(x,\) ip+(x,X) \ 



imply that there exists the limit 



^ ^ ^( X ,A) ^_(x,A) J ^ <p' a (x,X) j ^ A a (A) 
Furthermore, since expf— Qn\X)) — > I as n — > oo, there exists the limit 

M _ 1 / 1 1 \ / e"5 ar s^^'-)A Q (A) 

2 

It follows that |y4 a (A)| 2 = lim v an (e) /4, and substitution of this to the Weyl- 
Titchmarsh type formula (JHJ) completes the proof. □ 

The summary of this section is that the study of the spectral density of C a is reduceable to 
the study of the system (TT9]) and, more precisely, to the study of the behavior of lim v an (e) 

for small e. 

In the general case [y cr = Vj t ±) denote the coefficients (3, which may be different at 
different resonance points, as j3j t ±. Explicit calculations give: f3j t ± - ' ^" '="•''- '' 



2aW{V>+K±),<M^,±)} 

Remark 1. Coefficients (3j y ± are not necessarily non-zero, because they are proportional 
to the Fourier coefficients of p + (-,i / j t ±), which might be zero. E.g., consider the case of 
zero periodic potential, q(x) = 0. In this case one can choose the period a arbitrarily and 
the result is independent of the choice (except for the case u e ^j' Tor an U fixed a, the 

half-line is divided into spectral bands with coinciding endpoints, (^f) 2 ; (~7^) 2 , j > 0. 

The quasi-momentum is k(X) = aV~X, the critical points are = (- (j ' + 1 — {~})) 2 

and u jt - = (I (j + {f})) 2 , j > 0. Bloch solutions are iff±(x,\) = e ±ik ^~* = e ±Vlx and 
their periodic parts p + (x,X) = p~(x, A) = 1 have only one non-zero Fourier coefficient. 
Explicit calculation shows that = and /3j t - = for every j ^ \_^\ ■ Tor the single 

existing resonance point one has: vysuL^.. = u 2 and = M;- This coincides with the 

classical results on the Wigner-von Neumann potential [SO]- Our result concerning zeroes 
of the spectral density in this case is in accordance with the result of Hinton- Klaus- Shaw 



4. Model problem 

In this section our aim is to study the dependence on e, which can be arbitrarily small, 
of the limits of solutions to the system ({19]) . As we have shown in the previous section, 
this is equivalent to the study of the behavior of the spectral density of C a near critical 
points. 

Let us make few comments on the structure of the coefficient matrix of the model system. 

£n rr l f? ( t~\ n _ ( C0S V sinyp 



One can write it as follows: I + -D_ F „o~?, + FLfe), where D M = ( . is the 

« nw ' v \ — sin ip cos <p J 

matrix of rotation by the angle ip and 03 = ( J ^ j (notation for the Pauli matrix) 

is the reflection matrix. The presence of the latter plays a very important role. With 03, 
the system is elliptic for e G £^\{0} and hyperbolic for e = (we call the system elliptic, 
if it has a base of solutions of the same order of magnitude and hyperbolic in the opposite 
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case). Without 03, there is no change of type of the system at the point e = 0, the system 
is always elliptic. If 03 is absent, one can factor out the diagonal term of the first two 

summands 1 + j3 cos( - £n ^ and the coefficient matrix reduces to I + - sm( - £n ' > ( ® J -f Rn{e) 



v 1 / 

with some other uniformly summable sequence {i? n (e)}^L 1 . The matrix in the second term 
here is constant and diagonalizable. Its spectrum is purely imaginary, and the Levinson 
theorem immediately gives the answer (=asymptotics of solutions), which is uniform in 
£ G U including the point e = 0. We could say that the problem is "scalarized" in this 
case. The situation is different in our case, since the difference between the two eigenvalues 
is not pure imaginary: a (I + ^D_ en a 3 ) = jl ± -}. This leads to serious troubles in the 
analysis and exhibits a new phenomenon. 

One may expect that for sufficiently small values of e the magnitude of solutions is 
determined mostly by diagonal elements 1 ± (3 cos ^ £n ) . We want to transform the system in 
a way such that for every e G U (including e = 0) the limit of every solution as n — > 00 
exists (for the system (fl9|) this is not true: if e = 0, then one of the solutions grows as n 13 ). 
To this end, we make the following substitution: v n = exp {^[3 J™ cos ( £r ) c ; r j Un This leads 
to the system u n+ i = B n (e)u n with 



(22) B n {e) := exp 1-/3 / y -^dr 



n+l cosier) 



■ P(coB(en) sin(sn) 
n \ smieri) —cos{en) ' 



The existence of the limit of any solution of the system u n+ \ = B n (e)u n is equivalent to the 
convergence of the infinite product n^Li B n (e) (in fact for e 7^ this follows from Lemma 
H]). Moreover, most of the statements that we make about the asymptotic behavior of 
solutions of the system u n+ \ = B n (e)u n can be formulated in terms of products of matrices 
B n (e). In what follows we will choose the way of formulation depending on the convenience 
of its use. We are going to work with particular solutions determined by fixing their initial 
values, this is a discrete analogue of the Cauchy problem. To this end, we introduce the 
following notation: for given e G U and the vector of initial data / G C 2 define the vector 
sequence {u n (e, f)}™ = i by the recurrence relation 

u n+1 (e,f) := B n (e)u n (e, f), n > 1. 

Note that due to the decomposition of the exponential term the matrix B n (e) can be 
written as 

B n (e) = I+?-( ™(en)- f^cos(er)dr sin(sn) \ + ^ 

n \ sin(enj — cos{en) — J n cos{er)ar J 

with a sequence {R n (e)}^ =1 G l l (U). One can rewrite the system u n+ \ = B n (e)u n as 



(24) u n+ i - u r 



(3 ( cos(en) — cos(er)dr sm(en) 



sm(en) — cos(en) — J n+ cos(er)dr 



+ Rn(e) 



u. 



The behavior of the solutions can be observed in a scale of the variable y = n\e\ ("slow 
variable"). If one puts z(y,e,f) := u |x|( e > /) anc ^ divides by |e|, then ff24l) becomes 
approximately 



^ z{y + \e\) - z{y) 



P signe ( siny \ 

y \siny -2 cosy J \ £ \ LrJ 



z{y) 



The remainder |^|-R^jlJ ( £ ) is a step- wise constant matrix- valued function which is com- 
pressed in jiy times in the horizontal scale and stretched in ^ times in the vertical scale, 
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therefore it concentrates near the origin of the variable y and its L\ norm is preserved. So 
we may expect that the remainder term will be absorbed into a new boundary condition 
of the limit problem. Expression on the left hand side of ([25]) becomes the derivative as 
e — > 0, and one "obtains" the following equation for the limits h±(y, f) := lim z(y,e, f): 

e-s>±0 

(26) /4(v) = ±-( -° 7 y Wv), 

v 1 ±Kyj y y smy -2cos?/ J Kyj ' 

cf. (1540 . The remainder ^i?^ * j (e) plays a role only in determining the initial values 

h±(0) for the solutions of (1261) . It even suffices to know only R n (0) to determine this initial 
value, under the condition of continuity of R n {£) for every n. Thus the picture of the whole 
phenomenon can be described as follows: there exist two scales: "fast", discrete, n G N, 
and "slow", continuous, y G M+. The system first moves along the first ("fast") scale 
with e = 0. The limit of the solution as n — > oo for e = serves as initial value for the 
differential equation in the second ("slow") scale. Our aim in this section is to prove the 
following result which gives the exact formulation of the above considerations. We prefer 
to write an integral equation in the slow variable instead of the differential one, because the 
first has a unique solution while with the second one can have troubles due to the different 
behavior of the solutions near the origin in different cases depending on the value of /3. 

Theorem 1. Assume that functions R n {-) are continuous in U for every n G N, the 
matrices B n {e) are invertible for every n G N, e G U and the sequence {R n (e)}™ =1 G ^(U). 
Then for every y > and f G C 2 there exist two limits 

h±(yj) := lim u\ y I (e,f), 
e->±0 LRJ 

which satisfy the following integral equations: 

(27) h ± (y, f) = Jim MO, /) ± / ( exp (_ 2 /j* ^ ds ) J ) ^M*> /)*■ 
Moreover, the following four limits exist and are equal: 

(28) lim lim u n {e,f) = lim h ± {y,f). 

e-*±0 n-4oo y— >+oo 

Additionally, the linear map : / i— > lim u n (0,f) has rank one. 

n— >oo 

Remark 2. 1. It follows that the linear map f i— > lim lim u n (e,f) also has rank one. 
2. Note that lim u n (0, f) = lim lim u n (e,f) and lim lim u n (e,f) are the limits of the 

n— >oo n— »oo n— >oo 

same expression taken in the different order and that they do not coincide. Our aim is to 
prove that the second limit of these two exists. It will follow then that the spectral density 
of C a can have zeros at critical points. In fact we can rewrite the second expression as 
lim lim u\ y_ i (e, f). Here the convergence of the first limit is uniform in e (unlike the 

e— >±0 y~ >+oo L|s|J 

convergence of lim u n (e,f ) ), and this makes it possible to change the order of limits as 



in 

Proof. We divide the proof of Theorem [1] into four steps. 

Step I. A priori estimate and uniform convergence of the tail for the matrix 
product. We start with few technical results concerning the system u n+ \ = B n (e)u n . 
These include uniform boundedness of its solutions (in both variables, n and e) and the 
uniform with respect to e convergence in the slow variable y. 
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Lemma 5. Let B n {e) be given by (1221 and {R n (e)}^ =1 G l l {U). Then for every e G U the 
product 

oo 
n=l 

converges. If matrices B n (e) are invertible for every n G N and e G U , then for every 
non-zero e the product is invertible, while for zero e it is of rank one. 

Proof. This follows from the discrete Levinson theorem [5] . Using the expansion 

/ r n+1 cosr \ f3 f n+1 ( 1 \ 

(29) exp I —j3 / dr 1=1 / cos(er)dr + 01 — ] as n — > oo, 



\ J n '/ n Jn / 

we rewrite each matrix of the sequence B n (e) in the following form (at places, we drop 
the argument e in order to simplify the notation, hoping that this will not lead to any 
confusion) : 

(30) B n = I + V^+R%\ 
where {R { n\s)}™ =1 G l\U) and 

(31) ^ (1) (£) ■= - ( cos ^ ~ £ +1 cos(er)dr sin(£n) \ 

n \ sin(en) — cos^n) — cos(er)dr J 

One has: 

qos £ — 1 sin £ 

(32) / cos(er)dr = sin(en) H cos(en). 



For e 7^ the conditions of Theorem 3.1 from [5] are satisfied. The named theorem 
yields the existence of a base of solutions of the system u n+ \ = B n (£)u n , which have the 
asymptotics 

q ^ + o(l) and ( ^ J + o(l) as n — > oo. 

This is equivalent to the convergence of the product Y\^ = iB n {£) and its invertibility. 
For £ = the matrix is reduced to 

B n — { q -y _ 20 I + -Ri 5 ^ , 



n 



and the base of solutions changes accordingly: 





J J + o(l) and — p 



1 



oil] 



as n — y oo, 



which follows from the discrete Levinson theorem Theorem 2.2]. The existence of such 
a base of solutions is in its turn equivalent to the convergence of the product n^Li B n {fS) 
to a rank one matrix (since the second solution goes to zero asn-} oo). □ 

Remark 3. The rank one matrix = n^Li-^i(O) defines the linear map 
f ^ lim^oo u n (0,f). 

Asymptotics of solutions of the equation u n+ i = B n (£)u n given above demonstrate the 
change of the system type. For non-zero £ the system is elliptic (i.e., its solutions have the 
same rate of growth), while for e = the system is hyperbolic (i.e., there exists a base of 
solutions, which have uncomparable magnitudes; this yields the existence of a subordinate 
solution) . 

The following lemma states the uniform convergence of the product of matrices B n {£) 
in the slow scale. 
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Proof. The sequence given by f[31j) has the following property: 

4/3 



(33) 



k>n 



< 



n sin 



for every 116N and £ G C/\{0}. 



This easily seen using the equality (j52"|) . elementary estimates pf^l < 1, | cos | 1 1 < 1 and 
Lemma [2J It enables to define the sequence 

00 



k=n 



Then 



5n = / + gf -g( 2 j 1 + jR (f). 

Following the ideas of [UJ|5], we want to consider the Harris-Lutz transformation: 
(I - QSi)~ 1s «(/ - Q^)- If n > g with, say, 



(34) 



c 2 := 8/3 sup 



eeu\{o] I sin ll 

then it is easy to see that the estimate ||Q„ (e)|| < 5 holds yielding the invertibility of 
(/ — (5^4.1). For such values of n by a straightforward calculation one has: 

(35) BW = (/ - gS,)- 1 (/ + Qf) _ gWj (j _ g£>) + (/ _ Q^^R® (J - gW) 

= /+(/- g&r 1 (gil - gi 2) )gi 2) + (/ - q%) ~ x b® (/ - gi 2) ) . 



Using the trivial bounds ||V„ (2) || < 2||V r r J i; ||||C# s; ||, < 3\\Rn'\\, ||Kl 1J |l < 'f and the 

bound |igl 2) (£)|| < -r^TT, one has 

nan. 



II 11^(2) II || D (6)| 



(6) 



=:it 

>( 5 )|| IIT/Wll 2 /9 



£ ||k ( %)+^%)||< £ 



16/3 2 |e| 



3 E IW 5) ( £ )|| — ^ as ?/ — ^ -(-00 



~« 2/ sm 5 

n>^ n>^ I 2 1 

uniformly in e G t/\{0}. A rather rough argument repeating the scalar estimates yields: 



kl 



< 



exp E HK (2) (£) + 4 6) (£) 



hence the assertion of the lemma holds, if B n is replaced by 1 . Then, coming back to 
the product of matrices B n {e): 



(36) J] B ^ = ( J - & II B £\e) ( J - ^fl j+i( e 



fe> FT 



Due to the estimate fl3"3"|) . g, ^ , (e) — >■ as y — > +00 uniformly in e G [A{0}. Therefore 
the convergence to the identity matrix in (1361) is uniform. This completes the proof. □ 
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The next lemma completes Step I and proves the uniform boundedness of all partial 
products of matrices B n (e). This lemma is rather non-trivial and plays an important 
role in the rest of the proof. It is only due to the choice of the scaling factor in v n = 
exp \ B J™ cos ( £r ) d r j Un that these products are uniformly bounded. 

Lemma 7. Let {R n ( £ )}^Li £ l X {U)- Then there exists a constant c 3 such that for every 
e G U and every n < oo 



(37) 



k=l 



< c 3 . 



Proof. Using the decomposition of the exponent (12"9"1) again one can rewrite the sequence 
B n in the following form: 



1 



B n {e) 





n+l 



exp -2(3 / 



cos(er) 




+ V^{e) + R«Xe) 



with some {R ( J\e)}™ =l G l\U) and 



(38) V^(e 







n 



n+l 



cos(er)dr 



— | cos(en] 
n 



-1 
1 

n+l 

cos (er)dr 



(3 ( Q,os{en) sin(£:n) 
n \ sm(en) — cos(en) 



1 
-1 



. f s 

H — sinfen) 
n 



1 

1 



We are going to perform the variation of parameters in the discrete equation 

1 



u n +i 





n+l 



exp -2(3 f 



cosM dr \ I +VS 3) (e) + R^(e) 



u. 



considering it as a perturbation of the equation 



«n+l 

This leads to the following: 
1 





n+l 



exp -2(3 J 



cos(er) 



dr 



u, 



(39) u n (s,f) 







exp ( -2(3 fl 

n— 1 / 



n cos(er) 



dr 



f 







j V exp (" 2/3 ^ 



n cos(er) 



dr 



Using the Gronwall's lemma and a simple estimate 



(40) exp -2(3 



cos s 



ds < exp -2(3 



3tt 

2 COS s 



V£ 3) (e) + R ( ;\e))u k (ej). 



ds < 3 W , if < t < y < oo, 



one gets the following bound for the solution. 



n-l 



[u, 



{eJ)\\<^e W [^\K\e) + R k 7) {e) 



fc=i 
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In other terms, 
(41) 



k=l 



<3^exp l^\\vi%)+li%) 



k=l 



Furthermore, the identity (1551) leads to the following estimate: 

5/3 

\\Vn ( £ )|| — ~~^\ £ \ ^ or ever y ^ G N and e G U, 

which can be easily obtained with the help of the equality ( 1321) and explicit bounds 
I cos e-i | ^ 1 i sin e < M Thus, if n < j| (where c 2 is defined by (J5j}), then 



<- i i sine 

— 1 e — 2 



(42) 



fc=i 



k=i 



including the case of e = 0, n = oo. 

The estimate (|4ip grows to infinity with n, therefore the tail of the matrix product 



ought to be considered separately. If e ^ and n > ||, one obtains from the equality 



=I + Vf> +F£ } , see fl5BJ: 



(2) , d(6) 



n ^ 



<^p| ||^ (2) ( £ )+4 6) (^ 



\ 



/ 



Estimating the exponent in the same way as in the previous lemma one gets: 



(43) 



*=LftJ+* 



< 



16/3 2 




c 2 


sin|| 



fe=l 



Since both expressions (l4"2~j) and ( 145]) are bounded uniformly with respect to e, B n = 
(I - Q { ^ +1 )B^{I - Q^)' 1 and WQ^W < 1/2, the product ]JL[f2.\ + i is bounded 

for n > p uniformly with respect to e. Assertion of the lemma follows. □ 



Step II. Rewriting the system u n+ \ = B n (e)u n in the form of a Volterra integral 
equation in the slow scale. Consider the equation (1591) . 



, 1 

n-l / I 

m cos(er) ^ 



n-l / J Q \ 

+ £ ( exp (-2B ^ dr ) J K'^ + fl ^>) A 



which is equivalent to the system w„+i = B n (e)u n . On this step we rewrite it in an integral 
operator form. Fix an e G ?7\{0}. Put again n = A and divide the sum in (159]) into 

two sums, which contain terms and V^ 3 -*, respectively. Then write the second sum of 
the two as an integral in a new variable r putting k = [t\. Doing this we get for y > \s\, 
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since a piece-wise constant function appears in the second sum: 
1 

(44) U | A |(e)=| / JftJ cosher) , \ 1/ 



'W {£)= y0 exp [-2(5 ^Mrfr) 

L^J- 1 / 1 \ 

+ E ( exp(-2/3i + f J ^ r U^ 7) (^^ 

/ 1 \ 

+ / ' ( exp(-2/^]^r) J ^J^W^ 

Scaling the variable of the integration r = A we write this as a Volterra integral equation: 

(45) z(j/) =g(y)+ f K(y,t)z(t)dt. 

Jo 

Here we denote by z, g and K the following piecewise-constant (on intervals of length e) 
functions: 

*(y,e,f) ■= u [jl\ ( £ J), if y > l £ l> 

(46) ,(»,«,/):= ^ (_ VJ JjLftJ ajy.N j / 



LtStJ-i / 1 

+ £ 



J ( exp (-2/3 /i# J ^r)rfr) J ^M*,/). if 2/ > l £ 



\ , (5) 



:= | o exp (-2/3^ J+i) ^ I if ^ * < l £ l 

For all the other values of y (and t, if it is present) define these functions to be equal zero. 
After we have done this, we can successfully use standard operator methods. Before doing 
that let us observe the point- wise convergence of the kernel K and the "free term" g of the 
integral equation f )45|) . It is easy to see from the definition 



V^(e) = P [ n+1 cos(er)dr( - 1 



n 

that 



cos{sn) sm{en) 
1 ! ' n V sin(5n) — cos(en) 



1 



y 

and hence for y > t > 

(47) K(y, t, e) ± ( exp ^ J J ^ as s ^ ±0. 
Additionally, for every y > 

(48) ^ )£ )^(J (/ + fj4 7) (0)n fc (0,/)^ ase^O. 



ZEROES OF THE SPECTRAL DENSITY 



19 



This follows from the uniform bound given by Lemma [7| and from the property 
{Rn\s)}^ =1 G l x (U) by Lebesgue dominated convergence theorem. 

Remark 4. It is very important that the limit in (I48)) coincides with the lim u n (0, f). The 

n— ¥00 

latter can be obtained from the following variant of discrete Levinson theorem [HJ Lemma 
4.4, case (b)]. 

oo . 

Proposition 2 ([IS])- Suppose that V"T < 00 an( ^ that there exist an M such that for 

k=i M 

m 

every m > n the estimate ]~[ |A/| > -h holds. Let u be a solution of the system 



l=n+l 
U n +1 



A, 







1/X r 



Rr, 



U, 



If Y[ \Xi\ — °°; then both sides of the following equality exist and the equality holds: 
i=i 



lim — — ■ — 

nAi 



1 




k=i 



k „ 

i=i 



In the case of our system w n+1 = B n (e)u n , one can successfully apply this argument, 
putting u n = nV(0,/), X n = {^f and R n = Ri 7) (0). 

Step III. Convergence in the slow scale. Now we establish the convergence as e — > ±0 
of the solution u(e) in the slow scale. We remind that in our notation i £ ) = z (Vi £ )- 

Lemma 8. For every y > and f G C 2 there exist two limits 

h±(yj) := lim z(y,e,f), 

e—>-±0 

which satisfy the following integral equations: 

1 

exp (-2(3 J t y ^ds) 

Proof. For every y > and e ^ define the operator rC yo (e) in the Banach space 
^((O^C 2 ) by the rule 



(49) h ± {y,f)= lim Mn (0,/)± f ( 



(3sint 



h±{t,f)dt. 



(50) K ya (e) : u(y) m- / K(y,t,e)u(t)dt, y G (0;y ). 

Jo 

Denote K(y,t, ±0) := lim K(y,t,e) and analogously define two operators JC yo (±0). We 

consider only the case e — > +0 here. The second case can be treated in the same way. 

First let us see that the operator JC yo (e) converges to JC yo (+0) as e —> +0 in the space 
of bounded linear operators in Loo((0; y ), C 2 ). It is enough to show that 



(51) 



ry 

max / \\K(y,t,e) - K(y,t,+0)\\dt as e +0. 
y&[0;y ] J 



Fix a A > 0. Since both kernels are uniformly bounded in all variables, there exists a 
?/i (A) < y such that for every positive e G U 

ry 

max / \\K(y, t, e) - K(y, t, +0) \\dt < A. 

3/e[0;yi(A)] J 
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For the same reason 



(52) max !" \\K(y, t, e) - K(y, t, +0)\\dt 



3/S[yi(A);yo] Jq 



max 

yE[yi(A);y 



\\K(y,t,e) - K(y,t,+0)\\dt + O(e) as e -> +0. 



Using the estimates 1 ^ OS£ = 0(e) and 1 — = 0(e 2 ) as e — > 0, one can write 

«<*•*•<)= Lxp^/;^^^) oj gs '"[|}' J) + °( £ » as ^ +0 

for y and t such that 2e < y < yo and e < t < y — e, where 0(e) is uniform with respect 
to these y and t. We remind that 

1 \ /3 suit 

exp (-2/3 Sf^s) J ~ ^ ' 

The mapping (y; t) \— > exp (—2/3 J cos ^ ds j is uniformly continuous on the compact set 
yi(A) < y < ?/0) < t < y (notice that it has a discontinuity at the point y = 
t = 0), the function t i— )■ is uniformly continuous in the interval [0;?/o]- Therefore 
maxt 6[e . 2/ _ E] ||if(2/,t,e) - K(y,t, +0)|| = o(l) as e -»■ +0 uniformly in y e [3/1 (A); j/ ]. To- 
gether with (152]) this means that for sufficiently small positive values of e the maximum 



+0) 



max 

ye[yi(A);y ] 



[" \\K(y,t,e)-K(y,t,+0)\\dt< A, 
./o 



and hence the maximum in (15T1) has the same property. Since A is an arbitraty positive 
number, the convergence of operators JC yo (e) follows. 

Secondly, as mentioned in Remark HI for every y > one has: 

g(y,e) — >■ lim u n (0) as e — > 0. 

We face an obstacle here: this convergence is only point-wise and is not in the norm of 
£oo((0; yo), C 2 ). The same is true for the convergence of z(y, e), which we are going to prove. 
There is no sense in directly applying the inverse of the operator / — /C W) (e). The proper 
object to consider is the difference z(y, e) — g(y, e), which converges in Loo((0; y ), C 2 ). To 
see this let us rewrite the equation (J35J) in the following form: 

(53) z(e) - g(e) = (/ - JC vo {e))- 1 lC Vo {e)g(e) 

(since JC yo (e) is a Volterra integral operator, (/ — Ky^e))' 1 exists). It suffices to prove that 
K yo (+0)g(e) — > K yo (+0) lim^oo u n (0) in Loo((0; y ), C 2 ). By a direct estimate we have: 

max / K(y,t,+0) (g(t,e) - lim u n (0)) dt < /33 2/3 / \\g(t, e) - lim u n {0)\\dt. 

The right-hand side tends to zero as e — > +0 due to the point-wise convergence and the 
uniform boundedness of g{t,e) by Lebesgue dominated convergence theorem. From ( |53i) 
we obtain: 

z(e) - g(e) -+ (J - /C yo (+0))- 1 /C s/o (+0) lim u n {0) as e ->■ +0 in ^((0; y ), C 2 ) 

n— >oo 

From the point-wise convergence of g(y,e) it follows that for every y > 

z(y,e)^ ((/ + (/-/C s/() (+0))- 1 /C,o(+0)) hm u n (0)) (y) =: h+(y). 
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Consider the integral equation (145]) again. Taking the limit as e — > +0 turns it into the 
equation for h + , fj4H . This is possible due to the uniform boundedness of the solution 
z and of the kernel K by Lebesgue dominated convergence theorem. This completes the 
proof. □ 

Remark 5. The assertion of the last lemma includes the case when lim^oo w n (0, /) = 0. 
In this case h±(y, f) = 0. This happens only for one particular direction of the vector 
/eC 2 , because the rank of is one. 

Step IV. Convergence as y — > +oo and changing the order of limits. It remains to 
show that every solution of integral equations (J1HD has a limit at infinity and to see that 
it is possible to change the order of limits in lim lim u\ y i (e, f). 

Lemma 9. The following two limits exist: lim h±(y,f). These limits are non-zero, iff 

y— >+oo 

lim u n (0,/)^0. 

n— >oo 

Proof. Consider the integral equation for h + : 

h + (y) = lim u n (0, f) + J o ( exp { _ 2f3 °jy CJ ^ ds) lj^h + (t)dt. 

Let as write it as a differential equation using the following factorization property of the 
matrix in the integral: 



1 



1 



exp (-2/3 J t y ^ds) ) I exp (2/3 fi> ^ds 



x 



1 

exp (-2/3 J t °° ^ds) 



Multiplying from the left by |^ ^ ^2/3 J°° coss ds^j J an< ^ ^ a ^^ n S ^ ne derivative in y 



one has: 



if/ \ ( smy \ , / x 2/3 cos?/ / > r 

hUy) - - . ) h+(y) = — - „ , limu n 0,/ . 

+ Kyj y \ smy -2 cosy J +Kyj y \0 1 1 



By Lemma the vector lim u n (0, f) G C 2 is either proportional to ( ) or is equal to 
zero and therefore 

(54) K(v) = -( ° t y Wv)- 

v 7 + vyy y V smy -2 cosy y +vy/ 

If lim u n (0, /) 7^ 0, then h + (y) is not identically zero. Every non-zero solution of (|54|) has a 

n— >oo 

non-zero limit as y — » +oo, which follows for example from Harris-Lutz results [14, Theorem 
3.1, p.85] (one should take there A(y) = ( ° j),V(y) = W; J),i2(y) = 
0). For the second limit (of h_(y,f)) the proof is absolutely analogous. □ 
The following lemma is a combination of definitions and previous results: Lemmas El El 

E 

Lemma 10. The following two limits exist: lim lim u n (e, f) and equal to 

E— >±0 71— >O0 

lim h±(y,f), respectively. 
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Proof. By definition of h±, lim h±(y,f) = lim lim u\y_\(e, f). By Lemma [6j 

y— >+oo y— >+oo £— >±0 L|s|J 

u \ y-\{ e i f) converges uniformly with respect to e G £/\{0}. Therefore there exist two limits, 
lim lim u\ y_\(e, f) which coincide with the limits lim lim u\ JL \(e,f), respectively. 

e— >±0 y— >+oo L|s|J y— >+oo e— >±0 L\s\j 

Clearly lim u\ v i (e, f) = lim u n (e, f) for e 7^ 0, which implies that lim lim u n (e, f) = 
im«| A |(£,/) = lim 



lim lim u^ w j(e, /)= lim h±(y,f). This completes the proof. U 



We are ready now to finish the proof of Theorem [TJ Limits lim u\ y_\(e, f) exist by 

£^±0 L|s|J 

Lemma [HI the equality lim lim u n (e,f) = lim h±(y,f) has sense and holds true by 

e->±0n->oo y— >+oo 

Lemma [TTJ1 The map G has the rank one by Lemma (see also Remark [3]). □ 

5. Zeroes of the spectral density 

In this section, we prove the main result of the paper by applying Theorem [TJ 

Theorem 2. Let q\ e Li(R + ) ; cu |- ; and p' a (X) be the spectral density of the operator 
C a . Let the index j be greater or equal to zero. If the solution ip a (y, Uj i+ ) of ^ is not a 
subordinate one, then there exist two non-zero limits 

/4(A) 



lim 



2|e| 



Analogously, if \p a {y,v^S) is not subordinate, then there exist two non-zero limits 

/4(A) 



lim 



\\-v jt - 



a|W{V + (^ i + ),V_(^, + )} 







Proof. We use the notation of Section [3J v cr , (3 instead of Vj t ±, /3j t ±. This will yield both 
claims of the theorem. By ( 12"U|) . the solution ip a (y, A) is related to the solution {f Qjn (e)}^ =1 
of ( fl9|) . In turn, the latter is related by the scaling 

v n = exp (j3 J" cos ( £r ) dr \ u n to the solution {u n (e, v at i(e))}™ =1 of the system u n+ i = 

B n (e)u n . We now apply Theorem to the sequence [TJto {u n (e, v at i(e))}™ =1 . Non-subordinacy 
of (fairer) means that lim u n (0, w ai i(0)) ^ 0. Due to continuity of Vi a (-) in U, Theorem[TJ 

n— >oo 

implies the existence of the following two limits 

lim lim u n (e,v at i(e)) ^ 0. 

e— »±0 n— too ' 

Taking into account the scaling v n = exp (^(3 J™ co& ( £r ) d r j Un means that there exist 

two non-zero limits lim lim \e\^v a ^ n {e). This in turn yields the assertion of the theorem 

due to the Weyl-Titchmarsh type formula (|2"T|) and the property of the quasi-momentum 
that its derivative is positive inside spectral bands. □ 

6. Discussion 

As mentioned in Introduction, our result is to be compared to the one of Hinton-Klaus- 
Shaw [15]. The difference in the classes of operators considered may seem significant: in 
[15] there is no periodic background and the non-summable part of the potential is an 
infinite sum of Wigner-von Neumann type terms. Nevertheless, both problems can be 
written after a few suitable transformations in virtually the same form, see system ([6]) in 
our case and the formula prior to (2.2) in [T5]. We are also able to consider finite or even 
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infinite sum of Wigner-von Neumann type terms and to reduce the problem to the same 
form by analogous transformations. The methods that we use also have common traits 
like usage of the slow variable y = \e\n = \e\x/a or work with the limit integral Volterra 
equation like (T27|) which actually represents an equation of the same type as (2.75) in [T5] 
but is written in a different, more explicit form. However, the approach that we develop in 
the present paper to our mind has several essential distinctions compared to the approach 
of [IS]. 

1. We use a discretization of the differential system and work with the discrete model 
system of a "simple" form ffl~9]) . In addition to bringing certain technical (inessential) dif- 
ficulties it makes our method more universal and enables one to consider a wider class of 
problems. As an example, we can analyze by a direct application of Theorem [T] the struc- 
ture of zeroes of the spectral density of the discrete Schrodinger operator with Wigner-von 
Neumann potential. This is by far not the only possible example. Using the solution of the 
model problem one can also consider differential Schrodinger operators with point interac- 
tions at integer points with strengths of interaction forming a sequence of the Wigner-von 
Neumann form. 

2. We do not reduce the problem to a scalar one and use the analysis of vector equations 
instead. This makes the analysis simpler and more transparent. As a consequence of 
this, the form of the limit equation ( 127|) is explicit and can be guessed from heuristic 
considerations (see ( |24|) -( l26|) ). 

3. We avoid fine technically involved estimates of oscillatory integrals and in fact our 
proof is divided into two parts. First (Step I), we prove rather rough a priori estimates 
by dividing the positive half-line into subintervals M + = (0; o-] U (o! +oo) for sufficiently 
large c 2 . Second (Steps II-IV), we deal with the slow scale and prove the convergence using 
operator techniques. This is the main distinction of our method which makes it sufficiently 
simpler than the method of [T5]. 

We expect that using the technique developed in the present paper one can treat a more 
complicated case of the Wigner-von Neumann perturbation of non-Coulomb type with 
slowly decaying power part. We plan to address this problem in the future. 
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